Spontaneous Symmetry Breaking in Two Coupled Nanomechanical Electron Shuttles 
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We present spontaneous symmetry breaking in a nanoscale version of a setup prolific 
in classical mechanics: two coupled nanomechanical pendulums. The two pendulums 
are electron shuttles fabricated as nanopillars [l] and placed between two capacitor 
plates in a homogeneous electric field. Instead of being mechanically coupled through 
a spring they exchange electrons, i.e. they shuttle electrons from the source to the 
drain 'capacitor plate'. Nonzero DC current through this system by external AC 
excitation is caused via dynamical symmetry breaking. This symmetry-broken current 
appears at sub- and superharmonics of the fundamental mode of the coupled system. 

PACS numbers: 45.80.+r,85.85.+j,77.65.-j 



Spontaneous symmetry breaking is one of the funda- 
mental ideas in elementary particle physics when unified 
gauge field theories are applied [2[, e.g. breaking the 
gauge symmetry leads to the difference in the electromag- 
netic and the weak interaction. Commonly, one refers 
to a classical mechanics example to visualize symmetry 
breaking, such as a thin bar under pressure from both 
its clamping points. Obviously, the equations of motion 
are invariant around the symmetry axis along the beam. 
Once the force reaches a critical level the compressed 
bar will buckle and thus break the initial symmetry. Al- 
though, in the buckled state the bar will be in an en- 
ergetic minimum, the symmetry of this configuration is 
reduced. From a quantum mechanical point of view the 
key question is into which of the many degenerate en- 
ergetic minima the compressed bar will collapse to and 
what kind of fluctuation will have caused this transition. 
Naturally, for a large bar governed by classical mechan- 
ics it would be hard to associate the buckling with quan- 
tum mechanical states. The situation changes drastically 
when one considers quantum electro-mechanical systems 
(QEM) such as buckled nanomechanical beams which 
have been theoretically discussed by Carr et al In 
other words nanomechanical devices have the potential 
to address fundamental quantum mechanics [5|-[8j. 

In the following we will present electron transport mea- 
surements via two coupled electron shuttles at room tem- 
perature under AC/DC excitation. The two shuttles we 
fabricated have the same physical dimensions, are iso- 
lated from ground, and possess metallic tips. The actual 
realization of the device is shown in the scanning elec- 
tron microscope graph in Fig. 1(a): the two nanopillars 
forming the shuttles are placed in series between source 
and drain contacts. The diameters of the islands on top 
of each nanopillar is 65 nm with pillar heights of 250 nm. 
The inter-pillar distance is 17 nm and is much less than 
the gaps to source and drain, enhancing mode coupling. 



The distance from source to the first pillar and drain to 
the second one are identical, making the setup indeed 
very symmetric. Our particular focus will be on how 
breaking the local symmetry leads to the rectified current 
as recently conjectured by Ahn et al. [9j. The fundamen- 
tal outcome is that although both the nanomechanical 
and the nanoelectronic circuits are operated in the linear 
regimes an AC voltage can be rectified. Commonly, rec- 
tification only occurs based on a nonlinearity, however, 
such symmetry breaking in coupled nanomechanical pen- 
dulums is exactly what is predicted in Ref. 0|. 
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FIG. 1: (a) Aerial view in a scanning electron microscope of 
two coupled electron shuttles. The two nanopillars' gold tops 
have diameters of = 65 nm with a total nanopillar height 
of around 250 nm. The inset gives a top- view of the coupled 
nano-electromechanical pendulum. The scale bars are 200 nm 
in both graphs. Below: sketch of the measurement setup: an 
applied AC and DC bias voltage (Vbs) displaces the two pen- 
dulums by xi,X2 and the resulting current Ids is amplified, 
(b) Direct current Ids vs. bias voltage Vbs in the case of pure 
self-excitation, i.e. Vac = at room temperature. The upper 
inset shows low bias regime with the dashed lines indicating 
the occurrence of Coulomb charging effects (see text for de- 
tails) . The lower inset gives a finite-element simulation of the 
electric field distribution under an applied bias voltage (top 
view). The field intensity is maximal between the two pillars. 

The two nanopillars are milled out of a silicon-on- 
insulator (SOI) material, where the top crystalline silicon 
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is 190 nm thin and the insulating Si02 is about 350 nm 
thick. A 50 nm top gold layer serves as the electrical 
conduction path. The deposited metal serves as an etch 
mask in a dry etch step which mills out the SOI ma- 
terial around the pillars. We apply a CF4 plasma etch 
step and etch into the Si02 insulating layer, thus ensur- 
ing electron transport via the metallic islands (further 
processing details are given in Ref. [1Q|). The nanoscale 
circuit is placed in an impedance matched transmission 
line in order to minimize signal loss along the line. All 
measurements are performed under vacuum in a probe 
station at room temperature. The station is placed in a 
Faraday cage and is equipped with radio frequency con- 
tact probes covering the range from DC to 50 GHz ( 
a bias-tee allows AC/DC superposition with high preci- 
sion) . The equivalent circuit diagram is given in the lower 
part of Fig. 1(a): the two pillars can be individually dis- 
placed by x\ and #2, which lead to tunable resistances 
and capacitances (arrow boxes). 

In the first measurement we trace the direct current 
at different voltage bias values at room temperature as 
shown in Fig. 1(b). At first glance the /V-characteristic 
appears to be linear with only a small deviation (oc V 2 ) 
at above 1 V. Electron shuttling is leading to the non- 
zero current, which is caused by self-excitation (SE) as 
we have shown earlier for single nanopillars This 
can also be shown for lateral electron shuttles jll[l2j. In 
this particular case we assume that several mechanical 
modes of the coupled nanopillars (see below in Fig. 2) 
support the current. The ease of excitation is supported 
by the finite-element numerical simulation shown in the 
inset of Fig. 1(b). Obviously, the electric field intensity is 
dramatically enhanced around the two pillars, enabling 
the onset of mechanical oscillations under an applied DC 
bias voltage. Closer inspection of the low- volt age regime 
reveals the onset of Coulomb blockade (CB) in these seri- 
ally coupled electron shuttles, as predicted by Gorelik et 
al. [l3| . The dashed lines in the inset indicate the individ- 
ual Coulomb staircase, which is getting weaker towards 
voltages above 100 mV. A brief investigation of the trace 
reveals that the CB charging energy is of the order of the 
room temperature energy at 26 meV. Since we are op- 
erating the coupled shuttles at large AC-voltage swings 
(> 100 mV amplitude) as compared to the thermal back- 
ground (Vth — 26 mV), we can focus on the overall linear 
/^-characteristic. For all measurements presented here 
the gate voltage potential was left at ground. 

A radio frequency signal in the range of 1 MHz to 
1 GHz is applied via high-frequency probes in order to 
trace the mechanical spectrum of the coupled nanc-pil- 
lars, similar to earlier work on single nanopillars pi Hoj|. 
The resulting output current from the coupled electron 
shuttles is a time-averaged DC signal leading to the spec- 
trum shown in Fig. 2. The DC response of the coupled 
shuttles to the AC excitation is plotted at a DC bias 
voltage of Vbs = mV. A broad set of resonances can 
be identified representing different modes of the coupled 
nanopillars. On closer inspection of the overall spectrum 
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FIG. 2: (a) Spectrum of drain current Ids vs. / the ex- 
citation frequency of the coupled nanopillars. The current 
resonances are labeled according to the commensurate mode 
numbers p/q where the fundamental mode is found to be at 
cj /(27t) = 504 MHz from finite element simulations. The 
externally applied power of the radio frequency source was 
P — 4 dBm. Note the seven distinct current peaks in 
the reverse current (blue triangles). This pumping electro- 
mechanical pumping mechanism pushes electrons against the 
DC bias of Vbs = mV. The low-frequency region (★) is 
plotted in detail in (b). The red dashed lines are fits to the 
resonances taking into account different line shapes, (c) Fi- 
nite element simulation based on COMSOL showing the fun- 
damental mode of the coupled pendulums in a 3D (left) and 
top view (right). The color coding indicates the displacement 
within the nanopillars. 



we find that the multitude of peaks arises due to a para- 
metric instability caused by dynamical symmetry break- 
ing, as predicted by Ahn et al. [9]. The ratios noted at 
each peak follow (p/q)uo/(2ir) with p,q being integers. 
From finite element simulations we find that the most 
likely fundamental mode of the coupled nanopillars is at 
around uj /(2tt) = 504 MHz at p/q = 1, see Fig. 2(a). 
Unlike the theoretical prediction [9[, we find both even 
and odd numbers for q. The hierarchy of resonances is 
also slightly different from the theoretical result in that 
the current spectrum of the coupled electron shuttles has 
a richer structure with higher orders for p. We can at- 
tribute these differencies to more degrees of freedom in 
the real device allowing mechanical modes in all three 
spatial dimensions. 

With increasing external AC power we observe peak 
broadening, overlaping with neighboring peaks. We find 
the modes with the winding numbers, p/q generated by 
'Farey sums' of neighboring modes' fractions when over- 
lapping on the left and right: The numerator is the sum 
of the numerators of the fractions on either side, and the 
denominator is the sum of the denominators of the frac- 
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tions on either side [14[. For instance, p/q = 1/7 = (2 + 
2)/(15 + 13), 2/13 = (l + l)/(7 + 6), 2/11 = (l + l)/(6 + 5) 
and so on. A similar behavior was observed in nonlinear 
systems showing structure of Arnold tongues [15]. The 
appearance of fractions of the fundamental mode and 
Farey trees with neighboring modes' fractions indicate 
that we detect dynamical symmetry breaking. 

Most peaks of the spectrum show a forward current, as 
one expects under forward bias, while some of the peaks 
give a reverse current (blue triangles). This remarkable 
feature of the coupled shuttles, i.e. reverse shuttling or 
pumping, marks a bistability depending on initial condi- 
tions, such as positions and velocities of two islands at 
a given time. We need to note, as will be underlined 
theoretically below, that the resulting rectified current 
at zero bias observed here is solely due to a dynamical 
symmetry breaking and not by a mechanical or electronic 
nonlinearity. In detail the resonance shapes of the two- 
nanopillar system are analyzed in Fig. 2(b): a magnified 
part of the spectrum at below 100 MHz is shown with 
three fits (red lines) indicating the different line shapes. 
The resonance around 69 MHz follows the conventional 
Lorentzian, while the trace at 84 MHz shows a maximum 
and minimum. 
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FIG. 3: (a) Current response on- and off-resonance: at 
589 MHz (5) with an AC radio frequency excitation of -3 dBm 
a strong mechanical response of the coupled pendulums is ob- 
tained, while at off-resonance (77) only a weak current shows, 
(b) Bias dependence of the coupled nanopillars on-resonance S 
and off-resonance 77. While we find an ohmic response for both 
cases the red (5-trace is offset from the blue 77-trace. This off- 
set results from the rectification through the coupled nanome- 
chanical shuttles. 

We now turn to the details of the rectification effect 
due to symmetry breaking: as an example we select the 
pronounced resonance at (13/11) x ujq. The particular 
peak around 589 MHz is plotted in Fig. 3(a) and marked 
as 5. The off-resonance background is labeled with 77. In 
Fig. 3(b) the /^-traces for these two frequency points are 
given: as expected for the background the 77-line shows 
a linear dependence crossing through the origin of the 
graph. In contrast the £-trace reveals a clear offset un- 
der finite AC excitation at 589 MHz. This motion of 
the coupled charge shuttle indicates the existence of a 
bifurcation point [9]. 

The theoretical underpinning for this bifurcation of the 



coupled shuttle leads to the rectified current [9]. A si- 
nusoidal electromagnetic excitation Vq sin(co>7;) is applied 
to the source contact, inducing mechanical displacement 
of the shuttle at their eigenfrequencies. We can assume 
that the capacitances between the electrodes are con- 
stant, while the tunneling resistance R varies with elec- 
trode distance as Ri(xi) = i?i(0)e^ 1//A , ir^^i — x 2 ) = 
R 2 (0)e^- Xl ^ x ,R 3 (x 2 ) = R 3 (0)e- X2/x , where A is the 
tunneling length. The Xj refer to the three junctions 
(j = 1, 2, 3) in series, i.e. from source to the first nanopil- 
lar, in between the pillars, and finally from the second 
pillar towards drain. The coupled nanopillars can be as- 
sumed to show a linear response under modest excitation 
as measured in Fig. 1(b) leading to a mechanical displace- 
ment, which can be described by a differential equation 
in the adiabatic limit |9j 
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where x = x\ — x 2 is the relative coordinate, 7, uo 
denote the friction coefficient and the natural angular 
frequency of the shuttles, m gives the nanopillars mass, 
and L is the source-drain distance. There is a pair of 
bistable solutions, x(t) and —x(t). The time-averaged 
direct current is found to be 
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where we assume a symmetric configuration: Ri(0) = 
0.5^(0) = Rs(0) = R. For a single nanopillar with a 
symmetric source-drain the direct current with AC input 
voltage will always average out to zero, making the above 
integral zero. However, for two coupled nanopillars the 
symmetry will be broken leading to a DC signal as we 
observe in Fig. 3(b). We note that this current only can 
be found for two coupled nanopillars. 

In order to further verify the dynamical symme- 
try breaking and the resulting rectification we deter- 
mined the AC power dependence of the low-frequency 
end at around 80 MHz, revealing the expected Arnold 
tongues [16|. Fig. 4(a) shows the detailed power depen- 
dence of the output current with line plots with the power 
level increasing from -10 dBm to +4 dBm. In Fig. 4(b) 
the maximal current at 83 MHz is plotted vs. DC bias 
and AC power ramped up as the parameter. The in- 
set gives the scaling of generated direct current with re- 
spect to applied AC-power, indicating the rectification. 
In Fig. 4(c) finally the total power dependence is shown in 
a color scale plot: the boundaries of each tongue are seen. 
The experimental bifurcation diagram shows the grad- 
ual appearance of higher order resonances as the input 
power increases. The thin dashed lines trace the increase 
in width as expected for Arnold's tongues. The appear- 
ance of Arnold's tongues is also observed in mathmatical 
simulations of two charge shuttles coupled via electron 
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FIG. 4: Arnold tongues: (a) mechanical resonances at around 
80 MHz under increased AC power, (b) The DC bias depen- 
dence with gradually increased AC power, from -20 dBm (bot- 
tom) to -15 dBm (top). Inset shows the linear dependence of 
the rectified direct current upon increasing AC power, (c) 
Color scale plot of the total power dependence, revealing 
Arnold's tongues with the whited dashed lines being guides 
to the eye. Inset shows Devil's staircase plot depicting the 
width of each bistable region. 



tunneling as discussed in Ref. [9] . In each bistable region 
the fundamental mode and the external excitation are 
mode-locked in spite of the frequency mismatch. Devil's 
staircase depicts winding numbers of the mode-locked re- 
gions, as shown in the inset of Fig. 4(c). This staircase 
maps the widths of the Arnold's tongues, defined as the 
FWHM for a given power of AC excitation. Even though 
this two pillar system does not have a perfect left-right 
symmetry, driving the enhanced net current was possible 
due to the robustness of the observed mode locking. 

In summary we find dynamical symmetry breaking in 
coupled linear nanopillars traced in the current response 
of the pillars to an AC excitation. The identified reso- 
nances were agreed with the theoretical expectation by [9] 
except for some discrepanies due to the three-dimensional 
nature unlike the calculation model. The power depen- 
dance manifests the Arnold tongue structure where the 
bist ability of the system allows for nonzero currents. This 
setup has potential as a bifurcation amplifier for signal 
processing applications pj) and c an p rove to be useful 
for energy scavenging applications |18j . 
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